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Abstract 

The generalized Gerasimov-Drell-Hearn (GDH) sum rule is known to be very sensitive to QCD 
radiative and power corrections. We improve the previously developed QCD-inspired model for 
the Q 2 -dependence of the GDH sum rule. We take into account higher order radiative and higher 
twist power corrections extracted from precise Jefferson Lab data on the lowest moment of the 
spin-dependent proton structure function T^(Q 2 ) and on the Bjorken sum rule r^ _n (Q 2 ). By using 
the singularity-free analytic perturbation theory we demonstrate that the matching point between 
chiral-like positive-Q 2 expansion and QCD operator product l/Q 2 -expansion for the nucleon spin 
sum rules can be shifted down to rather low Q ~ Aqcd leading to a good description of recent 
proton, neutron, deuteron and Bjorken sum rule data at all accessible Q 2 . 
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I. INTRODUCTION 



The problem of the nucleon spin structure and the peculiarities of its underlying QCD 
description has attracted a lot of attention over the recent years [H, 0] . In particular, this is 
due to an enormous progress in experimental studies of the spin sum rules at low momentum 
transfer Q 2 , from the very accurate Jefferson Lab data on the lowest moment of the spin- 
dependent proton structure function Y\(Q 2 ) and on the Bjorken sum rule r^ _n (Q 2 ) in the 
range 0.05 < Q 2 < 3 GeV 2 [3J. This data provided a good testing ground for combining 
both the perturbative and non-perturbative QCD contributions. 

Theoretical description of the nucleon spin structure functions g p ' n at large Q 2 relies on 
the Operator Product Expansion, and at moderate Q 2 their sensitivity to the radiative and 
higher twist power corrections becomes significant [4]. Due to such a sensitivity the transi- 
tion to the entirely non-perturbative Q 2 region is rather cumbersome. This transition was 
earlier addressed in the QCD-motivated model [B| for the Q 2 -dependence of the generalized 
Gerasimov-Drell-Hearn (GDH) sum rule [6| making use of the relation to the Burkhardt- 
Cottingham sume rule |7j for the structure function g 2 , whose elastic contribution is the 
main source of a strong Q 2 -dependence, while the contribution of the transverse structure 
function, = gi + <?2, is smooth. The successful prediction of this model was the distinct 
crossover point of the proton data for T?(Q 2 ) at low Q 2 ~ 200 - 250 MeV 2 . Its subsequent 
modification including radiative and power QCD corrections, made the description far 
more accurate, which was required by the increased accuracy of the data. 

Now we enter a new level of increasing experimental accuracy, obtained in the recently 
published proton JLab data [3j. They lie above the model inputs at Q 2 > 1.5 GeV 2 (while 
displaying quite a similar shape) due to a noticeable sensitivity of pQCD part of Y p,n (Q 2 ) 
and to poorly known higher twist contributions ^4,6,.., as well as the axial singlet charge do- 

Our present goal is to improve the model for the generalized GDH sum rule for proton 
and neutron using the values of the power corrections /i^e,.. and singlet axial charge do, 
systematically extracted from the JLab data 0, [Io| and by performing a similar program 
of the smooth interpolation between large Q 2 and Q 2 = 0. As we will see we are able to 
achieve a rather good description of the data at all Q 2 values. 

The JLab data were obtained in the low Q 2 region and, therefore, a special attention is 
needed to the QCD coupling in this domain. While the 1/Q 2 term in the OPE works at rel- 
atively high scales Q 2 > 1 GeV 2 , higher-twist (HT) power corrections 1/Q 4 , 1/Q 6 , etc., start 
to play a significant role at lower scales, where the influence of the ghost singularities in the 
coefficient functions within the standard perturbation theory (PT) becomes more noticeable. 
It affects the results of extraction of the higher twists from the precise experimental data 
leading to unstable OPE series and huge error bars [3] • It seems natural that the weakening 
or elimination of the unphysical singularities of the QCD coupling would allow shifting the 
perturbative QCD (pQCD) frontier to a lower energy scale and to get more exact infor- 
mation about the nonperturbative part of the process described by the higher-twist series 

In this investigation, in order to avoid the influence of unphysical singularities at Q = 
Aqcd ~ 400 MeV, we deal with the ghost-free analytic perturbation theory (APT) [ll| (for 
a review on APT concepts and algorithms, see also Ref. [I2I]). which was recently proven to 
be an intriguing candidate for a quantitative description of light quarkonia spectra within 
the Bethe-Salpeter approach [l^], as well as in the recent higher-twist analysis of the deep 
inelastic scattering data on the F 2 structure function [3]. For completeness, we compare our 
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results obtained with conventional PT and APT couplings and, finally, discuss the related 
uncertainties and stability issues. 



II. FORMALISM 

A. OPE regime Q 2 > A.q CD 

To recall the basic ideas of the approach let us consider the lowest moments of spin- 
dependent proton and neutron structure functions g±' n defined as 



T p { n (Q 2 



dx gi' n (x, 



(2.1) 



From now on, it is understood that the elastic contribution at x — 1 is excluded from the 
moments, since it is the "inelastic" contribution which can be matched with GDH sum rule. 

At large Q 2 the moments T^ n (Q 2 ) are given by the OPE series in powers of 1/Q 2 with 
the expansion coefficients (see, e.g., Ref. }15|). In the limit Q 2 3> M 2 the moments are 
dominated by the leading twist contribution, ij% n (Q 2 ), which can be decomposed into flavor 
singlet and nonsinglet contributions: 



±a 3 + l -aAE NS (Q 2 ) + ^a E s {Q 2 ) 



i=2 



Q 



2i-2 ' 



(2.2) 



where E s and E NS are the singlet and nonsinglet Wilson coefficients, respectively, calculated 
as series in powers of a s 16). These coefficient functions for rif = 3 active flavors in the MS 
scheme are 



E NS (Q 2 ) 

Es(Q 2 ) 



71 



- 3.558 



71 
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20.215 
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(2.3) 
(2.4) 



The triplet and octet axial charges a 3 = gA = 1.267 ±0.004 [17] and a 8 = 0.585 ±0.025 [18 
respectively, are extracted from weak decay matrix elements. As for the singlet axial charge 
cio, it is convenient to work with its renormalization group (RG) invariant definition in the 
MS scheme a = ao(Q 2 = oo), in which all the Q 2 dependence is factorized into the definition 
of the Wilson coefficient Es(Q 2 ). For detailed discussion of the higher-loop stability of the 
coefficient functions and prescriptions used in actual calculations, see Ref. (lo| . 

We address both proton and neutron spin sum rules (SSRs), and the singlet and octet 
contributions are canceled out in their difference T p — resulting in the Bjorken sum rule 
H 



rr n (Q 2 ) = ^e ns (q 2 ) 

6 



E 



2i-2 



Q 



(2.5) 



The unphysical singularities at Q ~ ^qcd in the PT series for the coefficient functions 
E.gjQ 2 ) (12 .4p and E NS {Q 2 ) (12. 3 \ strongly affect the analysis of the spin sum rules at low 
Q 2 [l0|. Their influence becomes essential at Q < 1 GeV where the HT terms start to play 
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an important role. The "soft-frozen" a s models are free of such a problem, thus providing 
a more reliable tool of investigating the behavior of the spin sum rules in the low-energy 
domain. 

The moments of the structure functions are analytic functions in the complex Q 2 plane 
with a cut along the negative real demonstrated in Refs. |20j,|2l|]. On the other hand, 

the standard PT approach does not support these analytic properties. The APT method 



11] gives the possibility of combining the RG resummation with correct analytic properties 



of the QCD corrections. The consequence of requiring these properties to hold in the DIS 



description was studied previously in Refs. |22l . |23 



Let us recall that the expression for Y p { n (Q 2 ) in the framework of the analytic approach 
is completely similar to the one in the standard PT (12.21) : 



^i,apt(Q ) 



1 

12 



±a 3 + 



pAPT 



(Q 2 ) + 



e£ pt (Q' 



/'2 



Q 



2i-2 



(2.6) 



The corresponding NNLO APT modification of the singlet and nonsinglet coefficient func- 
tions is 



E^ T (Q 2 ) 



E 



J NS 
APT 



(Q 2 



1 -0.318 Af\Q 2 ) -0.361 AfiQ 2 ) 
1 -0.318 Af\Q 2 ) -0.111 AfiQ 2 ) 



(2.7) 



where A k is the analyticized k-th power of three-loop PT coupling in the Euclidean domain 
and defined as 



1 

7T ./„ 



Im([ai n) (-<7, rif)] k ) da 
o~+Q 2 ' 



n 



3. 



In the one-loop case, the APT Euclidean functions are simple enough [11 



(2.8) 
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1 dA 



k (3q dL 



(2.9) 



Analogous two- and three-loo p le vel expressions are more involved. However, according to 
the "effective log" approach [2H in the region Q < 5 GeV one may use simple one-loop 
expressions (12. 9p with the effective logarithm L* : 



A?l 3 (L) A?£ = A\l 3 (L*) , L* ~ 2 ln(Q/A^ 



yi eff — 



0.50 A (3) . (2.10) 



Thus, instead of the exact three-loop expressions for the APT functions, in Eq. ( 12. 7p one 

can use the one-loop expressions (12.91) with the effective A parameter A mo d = whose 
value is given by the last relation ( I2.10p . This model was successfully applied for higher- 
twist analysis of low-energy JLab data in Refs. 0, [Io|, and also in the T decay analysis in 
Ref. 0]. Note also that the APT couplings are stable with respect to different loop orders 
at low-energy scales Q 2 < 1 GeV 2 [12], contrary to the standard PT approach. 
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The APT functions Ak contain the {Q 2 )~ h power contributions, which effectively change 
the values of the //-terms, when going from the PT to the APT framework. In particular, 
by subtracting an extra (Q 2 ) 1 term induced by the APT series for the Bjorken sum rule 



Vin(g 2 /A^ } \J Q \Q 



.a) 2 



9a , J 1 \ , . A£ , „f 1 



where x ~ 0.43 and ~ 0.18 GeV is the effective one-loop Aqcd parameter, we get the 
relation between fJ^X PT coming into the APT expression ( 12. 6ft and the conventional fi^~ n 
from Eq. ( I2~2l) : 



/^(lCeY 2 ) KJpt + x^' 



M 2 M 2 



(2.11) 



Along with the conventional PT scheme, we will also apply the APT approach based 
on Eqs. ( 12. 6ft and ( 12.71) to construct the improved model for smooth continuation of per- 
turbative expressions for T^' n (Q 2 ) and its non-singlet combination r^ _n (Q 2 ) down to the 
non-perturbative region Q 2 — > 0. 



B. "Chiral" regime Q 2 < A 2 QCD 

For the purpose of a smooth continuation of T^' n (Q 2 ) to the non-perturbative region 
< Q 2 < Aq CD [H|, we consider firstly the Q 2 -evolution of the integral 

?M 2 > \ 1 2 /'' 

HQ 2 ) = ~^UQ 2 ) = -^rj o dx 9l (x,Q 2 ), (2.12) 

which is equivalent to the integral over all energies of the spin-dependent photon-nucleon 
cross-section, whose value at Q 2 = is defined by the GDH sum rule ^ 



h(P) = -&, (2.13) 

where /t^ is the nucleon anomalous magnetic moment. Then, the function h(Q 2 ) can be 
written difference 

h(Q 2 ) = I T (Q 2 )-I 2 (Q 2 ), (2.14) 

where 

?M 2 f 1 9M 2 f 1 

!t(Q 2 ) = -T^T / dxg T (x,Q 2 ), I 2 (Q 2 ) = — - / dxg 2 (x,Q 2 ). (2.15) 
V Jo V Jo 

The well-known Burkhardt-Cottingham (BC) sum rule [7| provides us with an exact 
expression for I 2 (Q 2 ), in terms of familiar electric Ge and magnetic Gm Sachs form factors 
as 

W) = \,G m (q Y G ^-^ Q2 \ (2.16) 
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where \i is the nucleon magnetic moment. As a consequence of the strong Q 2 behavior of 
the r.h.s. of Eq. (I2.16p . we get for large Q 2 



i 

2 



g 2 (x } Q 2 )dx\ Q2 ^ oo = 0, (2.17) 
so I2 is much smaller than Ii for large Q 2 . Now from the BC sum rule (I2.16p . it follows that 

J 2 (0) = ^ + ^ e (2.18) 
where e is the nucleon charge. Then the GDH value (12.131) is reproduced with 

M0) = ^f- (2-19) 

To summarize, from the above equalities ( I2.16p . ( 12 . 1 7[) and ( I2.19p . we can conclude that 
the BC and GDH sum rules together, lead to positivity of It{Q 2 ) for all Q 2 in the proton 
case and a vanishing difference between I T {Q 2 ) and I\{Q 2 ) for large Q 2 . Thus, I T (Q 2 ) is 
a smooth and monotonous function, and it is possible to obtain its smooth interpolation 
between large Q 2 and Q 2 = 0. 



III. IMPROVED MODEL FOR SMOOTH INTERPOLATION OF I T (Q 2 ) 

To improve the agreement between the model predictions and the experimental data, we 
consider the general asymptotic expression 



1 1, pert IV 



2M 2 



1 



1 



1 



— ( ±a 3 + -a 8 ) E NS (Q 2 ) + -a E S (Q 2 ) + ^ 
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i=2 



Q 2 *- 2 



(3.1) 



where the nonsinglet E NS and singlet Eg coefficient functions are defined in Eqs. (12.31) and 
( 12 .4p . respectively. Then, the perturbative expression for It, defined above the matching 
point Qq, is 



I T ,pert(Q 2 ) = 6(Q 2 - Ql) [h,pert(Q 2 ) + h{Q 2 )] 



(3.2) 



where Ii tPer t{Q 2 ) is calculated from Eq. ( 13. ip . while / 2 (Q 2 ) i n known from the BC sum rule 
( I2.16p . The smooth interpolation to the GDH value at Q 2 = ( I2.19P is difficult and cannot 
be performed analytically. Following the procedure developed in Ref. jHJ, we instead make 
use of the smooth extrapolation of the perturbative expression (13. 2p . to the nonperturbative 
domain Q 2 < defining the polynomial in positive powers of Q 2 as 



N 



lT,nonpert(Q 2 ) — ©(Qo — Q 2 ) 



1 d n I T . 



pert 



2\n 



(Q - Qi 



2\n 



(3.3) 



where N is the number of derivatives, which is a free parameter of the model, together with 
the matching point Q = Qo, which have to be chosen to satisfy 



'T, nonper 



li A e 



(3.4) 
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In practice, the easiest way to solve the problem is to fix the number of derivatives N and 
then to vary the Qo value until the relation H3.4[) is satisfied. It is interesting to note that 
taking N = 1 does not allow for such a solution. 

Such a procedure can be considered as a matching of the "twist-like" expansion in negative 
powers of Q 2 and the "chiral-like" expansion in positive powers of Q 2 |5j, which is similar 
to the matching of the expansions in direct and inverse coupling constants. 

Once we have obtained the parameters N and Qq, then the all-Q 2 expressions for the 
moments Y\' n can be restored from It(Q 2 ) defined by Eqs. (13.11) - (13.31) . by using Eqs. (12 . 12[) 
and (12.141) . As we will see below, this can be done within both the standard PT and 
singularity-free APT in the same way, leading to rather similar curves for (Revolution, 
except that in the APT case the matching point Qq playing a role of the "pQCD frontier" in 
this interpolation scheme is noticeably shifted down to lower Q 2 scales (see the next Section). 



IV. HIGHER-TWIST ANALYSIS 



A detailed higher twist analysis of the recent Jefferson Lab data [3[ on the lowest moments 
of the spin-dependent proton and neutron structure functions T^ n (Q 2 ) and r^ _n (Q 2 ) in the 
range 0.05 < Q 2 < 3 GeV 2 was performed in Refs. j§, In particular, including only three 



terms of the OPE expansion /i^.s i n Eq. (12. 2p . a satisfactory description of the data has 
been achieved down to Q 2 ~ 0.17 GeV 2 in conventional PT and down to Q 2 ~ 0.10 GeV 2 in 
the APT. 

The lower Q 2 involved, the higher twist contribution is needed to describe the data. As 
was shown in Ref. [lOj, there is some sensitivity of fitted values of /i 4 to the minimal scale 
Qmin variations; namely, it increases in magnitude when one incorporates into the fit the 
data points at lower energies. This property of the fit was treated as the slow (logarithmic) 
evolution /^(Q 2 ) (and a Q (Q 2 ) in the singlet case) with Q 2 which becomes more noticeable for 
broader fitting ranges in Q 2 , as discussed above. Indeed, fit results for /i 4 with taking into 
account the RG evolution with Q = 1 GeV, as a normalization point become more stable 
with respect to Q m in variations. 

However, there is still a problem how to treat the evolution of higher-twist terms f^6^,..(Q 2 ) 
which again may turn out to be important when one goes to lower Q 2 , since the fit becomes 
more sensitive to very small variations of /^6,8,.. with Q 2 . Since the evolution of the higher 



twists fie,8 is still theoretically unknown, they can be taken as free parameters [26]. This 
procedure leads to rather small xbof ~ 0.1 since more free parameters come into a fit at 
lower Q 2 . On the other hand, by including /i6,8,.. into the fit, one observes only a small 
change in /i 4 (lGeV 2 ) which demonstrates its stability down to lower Q 2 . Taking this 
into account, in order to reduce the number of free parameters, in the current work we apply 
another fitting procedure and determine, first, // 4 and a$ at higher scale Q — 1 GeV. Then 
we extract applying the known QCD evolution for /i 4 (Q 2 ) and clq{Q 2 ) and fixing them 
at 1 GeV from the previous fit. The number of free parameters does not grow in this case. 
The fitting domain is restricted from below by Q m %n defined by the condition x 2 '/ ' D.o.f < 1. 
The corresponding results are listed in Table HI Due to unknown evolution of /^(Q 2 ), which 
tends to be quite noticeable at lower Q 2 , we do not go below Q m in and do not take into 
account /i 8 -term here. 

The advantage of the APT analysis is the infrared and higher loop stability of the radiative 
corrections, as well as the stability w.r.t. Aq C d variations, leading to the stability and 
convergence of the higher twist series extracted from the data. Indeed, as we see from 
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Table HJ in the APT case the applicability of the perturbative expansion (12. 6p is somewhat 
shifted down to lower Q 2 , due to the absence of Landau singularities (see also Ref. [10| and 
references therein). 



V. ALL-Q 2 SPIN SUM RULES 

In the perturbative expression (13. ip we take into account the two-loop perturbative cor- 
rection in the singlet E$ and non-singlet E^s coefficient functions, as well as the twist-4,6 
contributions discussed in the previous section. To explore the infrared sensitivity of the 
model of the smooth continuation to Q = 0, we used two different sets of higher twist terms 
(with // 4 and /i 4 6 , respectively) and the corresponding singlet axial charge extracted from 
the data above a certain minimal scale Q min - 




FIG. 1: Proton spin sum rule function T\(Q 2 ) with respect to the combined set of J Lab and SLAC 
data. Results are shown with an account of the twist- ^ term (left panel) and the twist-4,6 terms 
(right panel). Corresponding perturbative parts are calculated in the framework of conventional PT 
(dotted lines) and APT (dashed lines). All-Q 2 model function obtained by the smooth interpolation 
of I^(Q 2 ) is also presented in PT (dash-dotted lines) and APT (solid lines). 



In Fig. [I] we present the proton spin sum rule function T^(Q 2 ) obtained by the smooth 
interpolation of the perturbative part lT, pe r P (Q 2 ) to the non-perturbative region Q 2 — > as 
described in Section III. Calculations taking into account one twist-4 term (left panel) and 
two twist-4,6 terms (right panel) listed above are performed within the conventional PT and 
APT. In Fig. [2] we show the Bjorken sum rule function Ti~ n (Q 2 ) calculated at any Q 2 in 
the similar way as T^Q 2 ). 

The all-Q 2 model functions T^Q 2 ) and T^~ n (Q 2 ) in both versions of the perturbation the- 
ory (dash-dotted and solid lines) are rather close to each other demonstrating the agreement 
between the singularity-free APT analysis at lower Q 2 and the usual PT one at relatively 
higher Q 2 . Also, as one can see from the comparison of the left and right panels the results 
of the interpolation do not strongly depend on the number of higher twists included and, 
hence, on the border between perturbative and non-perturbative regimes. This exhibits 
a sort of duality between them implying that the experimental data in the wide intermediate 
region Aq CD < Q 2 ~ 1 GeV 2 can be described equally well either by OPE 1 / Q 2 -series or by 
"chiral-like" Q 2 -series. 

We studied the sensitivity of above results w.r.t. variations of the number of derivatives 
N in Eq. (I3.3P being the number of positive ~ Q 2% power terms. As mentioned above, at 
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TABLE I: Combined fit results of JLab and SLAC data on the Bjorken SR and proton SSR for the singlet 
axial charge arj, and the higher-twist terms \i± and \i§ defined at the normalization point Q 2 = 1 GeV 2 
in the APT and the standard PT approaches, along with the matching value Qq. Corresponding curves 
for F^Q 2 ) and r^~ n (Q 2 ) are shown in Figs. [T] and [21 respectively. Typical values of x 2 /D.o.f are close 
to unity. 



Method 


Target 


QLn, GeV 2 


a 


/x 4 /M 2 


«s/M 4 


Ql GeV 2 




p-n 


1.0 




-0.060(3) 





1.1(2) 


NLO PT 




0.3 




-0.060 


0.010(2) 


0.8(2) 




proton 


1.0 


0.34(3) 


-0.056(3) 





1.0(1) 






0.5 


0.34 


-0.056 


0.010(2) 


0.8(2) 




p-n 


1.0 




-0.058(3) 





1.0(1) 


NLO APT 




0.2 




-0.058 


0.010(1) 


0.3(1) 




proton 


1.0 


0.37(2) 


-0.063(2) 





1.0(1) 






0.3 


0.37 


-0.063 


0.011(1) 


0.3(1) 




FIG. 2: Bjorken sum rule function n (Q 2 ) with respect to the combined set of JLab and SLAC 
data. The meaning of curves here is the same as in Fig. d 

lower Q 2 we need more higher 1 / Q 2 -power twist terms. In the same way, going up from very 
low Q 2 we observe analogously that to describe the data at higher Q 2 we need more ~ Q 21 
power terms, i.e. a higher value N. 

The minimal number of derivatives N min , which is necessary to perform the smooth 
extrapolation according to Eq. (13. 3p in the conventional PT case and with one \i± term only, 
is N m i n — 4. Corresponding matching value between perturbative and non-perturbative 
domains in this case is found to be Ql = 1.0 ± 0.1 GeV 2 for the proton SSR and Q 2 = 
1.1 ± 0.2 GeV 2 for the Bjorken SR (see Table [I]). However, if one increases the number 
of Q 2 -power term up to = 6, the applicability of the "chiral-like" expansion raises up 
to Ql ~ 1.4 GeV 2 for the proton SSR and Q 2 ~ 1.5 GeV 2 for the Bjorken SR. Similar 
observation was made earlier in Ref. 

In the framework of APT the minimal number of derivatives N min = 3 is even smaller 
than in the conventional PT. In this case, if only one \x± term is included then the matching 
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value turned out to be same as in PT: Q$ ~ 1.0 GeV 2 for both the proton SSR and p — n 
demonstrating the similarity of the APT and PT predictions at Q > 1 GeV. 

However, analysis at lower Q 2 including two /i^ terms leads to quite different Qo values 
for PT and APT. In this case, for the proton SSR and Bjorken SR we have Qq ~ 0.8 ± 
0.2 GeV 2 (PT, N min = 3) and Q 2 ~ 0.32 ± 0.1 GeV 2 (APT, N min = 2). Such a shift 
of the border between perturbative and non-perturbative domains in the APT is a direct 
consequence of the disappearance of the unphysical singularities in the radiative corrections, 
and confirms the similar conclusion made in Refs. jl, 10 . 



-0.02 



-0.04 



-0.06 



-0.08 



-0.10 




Q„=0.3 GeV'(^PT) | 



Q„=0.8 GeV (PT) 



Q =0.S GeV (PT) 



Q =0.3 GeV (APT) 



FIG. 3: Neutron (left) and deuteron (right) spin sum rule functions, T™(Q 2 ) and T d (Q 2 ), with 
respect to the combined set of J Lab and SLAC data. Results are shown with an account of twist-4,6 
terms. The meaning of curves here is the same as in Fig. [71 



Finally, in Fig. [3] we show the neutron spin sum rule function T'KQ 2 ), which is simply 
obtained from the difference rf(Q 2 ) — T{~ n (Q 2 ), and the deuteron spin sum rule Tf(Q 2 ). 
We also present its perturbative PT and APT parts together with less precise data. Both 
versions of the perturbation theory predict monotonous curves for r"' d (Q 2 ) at any Q 2 . Com- 
parison between them and the results of Ref. 1(| demonstrates the consistence of the direct 
fits to the data and the predictions of the generalized GDH sum rule. 



VI. CONCLUSION 

In the current paper we have considered the all-Q 2 model for the generalized GDH sum 
rule, constructed by the smooth interpolation of It{Q 2 ) between large Q 2 and Q 2 = 0, in the 
framework of both the conventional PT and the ghost-free Analytic Perturbation Theory. 
We used the values of the power corrections /i4,6,.. and singlet axial charge do, systematically 
extracted from the precise JLab data. We achieve a rather good description of the proton 
data on T^(Q 2 ) at any Q 2 values. We also present an improved description of the neutron 
data, as well as the Bjorken sum rule data at all experimentally accessed Q 2 . 

The results of the smooth interpolation Tl(Q 2 ) and Tl~ n (Q 2 ) do not strongly depend 
on the number of higher-twist terms, and on the border Qq between perturbative and non- 
perturbative regimes. This exhibits a sort of duality between them implying that the exper- 
imental data in the wide intermediate region Aq CD < Q 2 ~ 1 GeV 2 can be described equally 
well either by OPE l/Q 2 -series or by non-perturbative "chiral-like" Q 2 -series. Within the 
analytic PT the "pQCD frontier" being the matching value between Q 2 - and l/Q 2 -power 
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series naturally decreases from 1.0 GeV 2 with single ^4 down to 0.3 GeV 2 with extra iiq- 
term included, which is significantly lower than the corresponding value in conventional PT 
Qq ~ 0.8 GeV 2 . Such a shift of the border between perturbative and non-perturbative do- 
mains in the APT is a direct consequence of the disappearance of the unphysical singularities 
in the radiative corrections. 
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